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ABSTRACT 

Recent work suggests blue ellipticals form in mergers and migrate quickly from the 
blue cloud of star-forming galaxies to the red sequence of passively evolving galaxies, 
perhaps as a result of black hole feedback. Such rapid reddening of stellar populations 
implies that large gas reservoirs in the pre-merger star-forming pair must be depleted 
on short time scales. Here we present pilot observations of atomic hydrogen gas in four 
blue early-type galaxies that reveal increasing spatial offsets between the gas reser¬ 
voirs and the stellar components of the galaxies, with advancing post-starburst age. 
Emission line spectra show associated nuclear activity in two of the merged galax¬ 
ies, and in one case radio lobes aligned with the displaced gas reservoir. These early 
results suggest that a kinetic process (possibly feedback from black hole activity) is 
driving the quick truncation of star formation in these systems, rather than a simple 
exhaustion of gas supply. 

Key words: galaxies: elliptical and lenticular, cD, galaxies: evolution, galaxies: for¬ 
mation 


1 INTRODUCTION 

Star-forming galaxies show a correlation between stellar 
mass and star formation rate whose normalisation varies 
with red shift, but whose scatter remains tight out to high 
redshift llNoeske et al.l 120071 : iPeng et ahl l2010l : lElbaz et al.l 
l201ll) . This correlation can be explained as a balance be¬ 
tween gas inflows fr om cosmological scales and outflow s 
driven by supernovae dBouche et al.ll2Q10l : iLillv et ahll2013l) . 
Quenching of star formation causes galaxies to depart from 
this steady-state along pathways that depend on the evolu¬ 
tion of the gas supply and reservoir (tSchawinski et aPl2014l) : 
if the cosmological gas supply to a galaxy is shut off, it will 


* E-mail:ivy.wong(9uwa.edu.au 


slowly exhaust its gas reservoir by forming stars at a declin¬ 
ing rate. If instead the gas reservoir fuelling star formation 
is destroyed effectively instantaneously, then star formation 
will cease and the galaxy will redden rapidly within 1 Gyr, 
as is observed in early-t y pe galaxies (ETG s; Kavirai et al.l 
l201ll : IWong et al.l 120121 : ISchawinski et all I20l4) . Interest¬ 
ingly, this timescale of 1 Gyr is consistent with that pro¬ 
posed for the transition of a merger-dri ven Ultr alu minou s 
Infrared Galaxy (ULIRG) to an ETG (lEmonts et al.ll2006l ). 
Also, it has been kn own for a while tha t ULIRGs have a 
high merger fraction (ISanders et al.|[l988l ). 


In addition, negative feedback from a galaxy’s central 
active galactic nuclei has often been invoked in galaxy for¬ 
mation models to slow do wn t h e star format ion history 
of simulated galaxies (e.g. ICroton et ali [iooiih . However, 
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Figure 1. Top row: optical gri SDSS llAhn et al.ll20lj) composite images of the four blue early-type galaxies in our sample. The galaxies 
are arranged from left to right in terms of NUV — u colour from blue to red as a tracer of time. The white scale bars represent 10 
arcseconds. Bottom row: zoomed-out optical images with the WSRT radio observations overlaid. The white contours (asinh-stretch) 
represent the Hi and the green contours (logarithmic stretch) the radio continuum. The lowest radio and Hi contours begin at 3er and 
the highest contour is matched to the peak flux density of 0.87 mjy beam " 1 . 1.51 mjy beam -1 , 0.86 mjy beam -1 , 44.79 mjy beam -1 
for J1237, J1117, J0900 and J0836; respectively. It should be noted that the nuclear 1.4 GHz continuum point sources are not shown for 
J0900 and J1117 to avoid confusion with the Hi contours. 


oDservational evidence tor sucn teed P ack is curren tly 
available for a few i n dividual galaxies llHarrison et al.l 
NvlandetalJ [2013]; |AlatalojtalJ 201ll: iHota et aiT 


Croston et al.l 20081 : Kharb et al.l 120061 ') . With the advent 


of very large surveys; much progress has been made to¬ 
wards understanding of the co-evolution between galaxies 
and their central AGN, in particular, the connection be- 
tween classical bulges and central supermassive black holes 
llHeckman fe Best1l20 1 4l : iKormendv fc Holl2013h . 


Blue ETGs do not fit the canonical bimodal scheme 
of red ellipticals versus blue s pirals and are unlikely to be 
the descendants of blue spirals dToieiro et al.|[2013l) . Rather, 
the blue ETGs appear to be transition-type galaxies and 
are probable predecessors of local post- starburst galaxies 
llSchawinski et ahll2009l : IWong et al.ll2012l ). Previous CO ob¬ 
servations of blue ETGs show that the molecular gas reser¬ 
voirs are being rapidly destr oyed during the compos ite star 
formation and AGN phase (ISchawinski et al.l 120091) . Local 
post-starburst galaxies are defined to be galaxies which have 
recently ceased forming stars. These galaxies typically show 
little Ha or [Oil] emission (indicative of current star for¬ 
mation), but have strong absorption line signatures indica¬ 
tive of a young stellar population. In addition, local post- 
starbursts consists mostly of galaxies which have similar 
early-type structural properties to red sequence galaxies of 
comparable masses (IWong et al.ll2012l ). 


Earlier attempts to map the Hi content of post- 
starburst galaxies only found Hi in one of five targets which 


consi sted of an interacting pair of galaxies (IChang et alJ 
l200ll) . All subsequent campaigns to detect Hi in post- 
starbursts have been optimised towards detecting diffused 
low-surface brightness gas, while compromising on the an- 
gular resolu tio n re quired to map the location of the gas (e.g. 
IZwaan et al.l[2013h . 


To determine the physical process responsible for shut¬ 
ting down star formation, we posit that post-starburst galax¬ 
ies are too evolved and that the smoking gun for quenching is 
more likely to be found in its predecessor population, namely, 
the blue early-type galaxies. 


Using the Westerbork Syntheses Radio Telescope, we 
study the Hi content of a pilot sample of four blue early-type 
galaxies—the progenitors to post-starburst galaxies because 
1) these galaxies are still currently star-forming and there¬ 
fore there should be gas where there are stars being formed; 
and 2) the state of the gas morphology and dynamics will 
shed light on why these galaxies will soon stop forming stars. 


Section 2 describes the blue early-type sample. The pi¬ 
lot observations and data processing methods are described 
in Section 3. We discuss the radio continuum and Hi imag¬ 
ing results in Section 4 and Section 5, respectively. Section 
6 provides a summary of our results. The AB magnitude 
system is used throughout this work. 
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Figure 2. Colour-stellar mass diagram of ETGs found from SDSS 
and Galaxy Zoo. The contours represent the number density of 
ETGs occupying the colour-stellar mass region known as the ‘red 
sequence’. All four pilot sample galaxies are located in the op¬ 
tical ‘green valley’ of the colour-mass diagram (as demarcated 
between the solid green lines). All galaxies show the presence of 
intermediate age stellar populations by their ‘green’ u — r colour. 


2 BLUE EARLY-TYPE GALAXIES 


We obtain the photometric and spectroscopic data from 
the S loan Digital Sky Surv ey (SDSS) DR7 l|York et al.1 


120001 : lAbazaiian et al. 


2009) for all objects classified as 
20021 ) . The sample of 204 low-redshift 


‘galaxy’ IlStrauss et al.1 
(0.02 < z < 0.05) blue ET Gs were identified using 


the Sloan Digital Sky Survey jAdelman-McCarthv et al 


200S, SPSS; ) and the Galaxy Zoo project (ILintott et al 


2008 ; ISchawinski et all 120091 : ILintott et, al.ll201ll ). The blue 


ETGs account for 5.7% of all ETGs found within the 
same redshift range and a re the most actively star¬ 
forming population of ETGs llSchawinski et ahl 12009). Fur¬ 
ther description of th e blue early-typ e gala xy selec¬ 
tion can be found in ISchawinski et alJ Q2009J) and at 
http://data.galaxyzoo.org. 


2.1 The pilot sample 

To investigate the fate of the gas reservoir when star for¬ 
mation is quenched, we select four ETGs with compa¬ 
rable stellar masses and UV/optical colours indicative of 
rapid quenching: all four lie in the optical u — r green val¬ 
ley, indicating they have significant intermediate-age stel¬ 
lar populations, and their NUV — u UV/optical colours 
range from very blue, consistent with ongoing star forma¬ 
tion (J1237+39), to the redder colours of passive, quenched 
galaxies (J0836+30)—see Figure [5] Due to the diffused na¬ 
ture of Hi gas, we have also selected the nearest galaxies 
which can be observed by the Westerbork Synthesis Radio 
Telescope (WSRT). The optical/stellar morphologies of this 
sample are presented in the top row of Figure [l] 

Stellar population modelling find two distinct star for¬ 
mation truncation timescales for early- and late-type galax¬ 
ies whereby the early-type galaxies appear to have fast 
quenching timescales (^1 Gyr), while the late-type galaxies 


Figure 3. The UV/optical colour-colour diagrams of blue early- 
type galaxies (green density contours) and the entire early- 
type galaxy population (black shaded contours). The model stel- 
lar population tracks with varying quenching timescales from 
ISchawinski et all d2014h are also overlaid. The large and small cir¬ 
cles along the evolutionary tracks represent 1 Gyr and 100 Myr, 
respectively. Two out of four pilot sample galaxies show an ab¬ 
sence of very young stellar populations indicated by red NUV — u 
colours. Those three galaxies also show either a displaced HI reser¬ 
voir (J1117-1-51) or a completely ejected HI reservoir (J0836+30 
and J0900+46). The galaxy which is bluest in NUV — u features 
a rotating HI disk. 


have quenching timescales of several Gyrs (ISchawinski et ahl 
120141 ). Figure [3] illustrates the different evolutionary stages 
represented by each of the four ETGs selected for this pilot 
study on a UV/optical colour-colour diagram. The colour 
of model stellar populations along different star formation 
quenching tracks are shown as solid lines where each large 
circle marks 1 Gyr and intervals of 100 Myr being shown 
by small circles along the stellar evolution tracks. The black 
shaded contours represent the low-redshift early-type pop¬ 
ulation and the green contour represent s the population 
of blu e early-type galaxies selected from ISchawinski et all 
(2009). The observed stellar populations of J0836+30 and 
J0900+46 appear to be relatively evolved along the fastest 
quenching evolutionary pathwa y that occurs on timesca les of 
several hundred Myr to 1 Gyr (ISchawinski et aPl2014l ). On 
the other hand, the stellar population colours of J1117+51 
and J1237+39 are less evolved and consistent with earlier 
stages of quenching. We note that J1117+51 and J1237+39 
appear to favour slower quenching pathways in Figure [3] 
However, this could be due to differences between the start 
times for quenching in these galaxies relative to the fiducial 
starting point for quenching in the models. The Hi proper¬ 
ties of each sampled galaxy appears to be consistent with its 
respective evolutionary stage. See Section 5 for more details 
of the observed Hi properties of this sample. 

I 11 addition, all four galaxies show optical emission lines 
dominated by activity other than star formation: J0836 and 
J0900 show Seyfert-like lines while J1117 and J1237 show 
low-ionization lines which could b e due to black hole ac¬ 
cretion, shocks or evolved stars dSarzi et al. I l2010l ). Figure [4] 
presents the Baldwin, Phillips & Televich (BPT) diagnostic 
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Figure 4. Optical nebular emission line diagrams. The nebular emission line ratios of the entire Schawinski et al. (2009) sample of 
blue early-type galaxies are plotted as black points. The blue solid points represent the objects selected in this pilot sample. The solid 
and dashed line are two lines used to separate the nebular emission originating from star formation rather than that of Seyfert or 
LINER activity (jKewlev et al.l[200ll: iKa uffmann et al.ll 2003lb The dotted-dashed line differentiates between ionisation levels obtained 
from LINER versus Seyfert cores ISchawinsk i et al. [ booL 


diagram (iBaldwin et al.lll98ll l which shows the emission line 
ratios observed for blue early-type galaxies. Table Q] lists the 
optical properties of this pilot sample. Stellar masses were 
obtained from spectral fitting (ISchawinski et al.ll2009l 'l. 


3 PILOT SURVEY 
3.1 Radio observations 

We performed imaging of the 1.4 GHz radio continuum and 
the atomic Hydrogen (Hi) emission of our pilot sample using 
the WSRT in the Netherlands between June and November 
2012. As the WSRT is an east-west interferometer, our ob¬ 
servations were divided into several epochs to optimise the 
Mv-coverage. The single 20 MHz band is sampled over 1024 
channels in the maxi-short baseline configuration in order to 
obtain a spectral resolution of 4 km s -1 spanning the range 
of 4000 km s -1 . 

Using this setup, a 24-hour on-source integration should 
result in an RMS noise level of 0.54 mjy beam -1 at full an¬ 
gular resolution over a FWHM of 8.25 km s -1 in velocity— 
corresponding to a column density sensitivity of 1.46 x 10 19 
atoms cm -2 per resolution assuming a uniform weighting]. 
Due to timetabling constraints, we were not able to obtain 
24-hour on-source integration for each of the four galaxies. 
Total integration times for each galaxy, the resultant synthe¬ 
sised beam properties and sensitivities are listed in Tabled 


1 Noise estimations are obtained from the 
WSRT exposure time calculator found at 
http://www.astron.nl/ oosterlo/expCalc.html 



Figure 5. Multicolour composite of the J0836+30 field (centred 
on J0836+30) where the background colour image is a multicolour 
composite of the five SDSS optical ugriz bands. The 1.4 GHz 
radio continuum and Hi emission from these pilot observations 
are shown in magenta and cyan, respectively. The yellow ring in 
the bottom-left of the figure represents the beam of the radio 
observations. 


3.2 Data processing 

The observations are calibrated and processed in the 
standard manner using the Miriad Reduction software 
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Table 1. Optical properties of our low-redshift blue early-type galaxy sample. 


SDSS ID 
(1) 

Galaxy 

(2) 

RA (J2000) 
(3) 

Declination (J2000) 
(4) 

Redshift 

(5) 

Distance 

(6) 

r 

(7) 

u — r 

(8) 

log(AU) 

(9) 

587735044686348347 

J0836+30 

08:36:01.5 

+30:15:59.1 

0.02561 

105 

-21.1 

2.27 

10.7 

587731887888990283 

J0900+46 

09:00:36.1 

+46:41:11.4 

0.02748 

113 

-21.3 

2.20 

10.5 

587732135382089788 

J1117+51 

11:17:33.3 

+51:16:17.7 

0.02767 

115 

-21.4 

2.33 

10.6 

587738946685304841 

J1237+39 

12:37:15.7 

+39:28:59.3 

0.02035 

84 

-20.9 

2.16 

10.3 


Col. (1): SDSS object identification. Col. (2): Galaxy identification used in this paper. Col. (3): Galaxy center’s right ascension. Col. 
(4): Galaxy center’s declination. Col. (5): Redshift. Col. (6): Distance in Mpc. Col. (7): SDSS r-band magnitude. Col. (8): Optical u — r 
colour. Col. (9): log (stellar mass) in solar mass as estimated by spectral fitting. 


Table 2. WSRT observation summary. WSRT observation integration times and resulting beam properties are detailed for each of our 
pilot blue galaxies. 


Galaxy 

(i) 

Observations dates 
(2) 

Total integration time 

(3) 

Beamsize 

(4) 

Beam Position angle 
(5) 

RMS Hi 

(6) 

RMS 1.4 

(7) 

J0836+30 

Nov 15,25 

20.3 

29.5, 13.9 

-2.0 

4.1 

0.083 

J0900+46 

Oct 23, Nov 11 

10.3 

24.0, 12.1 

-1.7 

11.0 

0.062 

J1117+51 

June 30, Nov 18 

24.0 

19.8, 14.5 

9.6 

3.7 

0.082 

J1237+39 

Aug 24, Oct 25, Nov 21 

24.2 

27.6, 17.5 

4.7 

32.6 

0.038 


Col. (1): Galaxy identification. Col. (2): Observation dates. Col. (3): Total on-source integration time in hours. Col. (4): Synthesised 
beam diameter in arcseconds. Col. (5): Synthesised beam position angle in degrees. Col. (6): RMS noise level for the Hi emission in 
mjy beam -1 . Col. (7): RMS noise level for the 1.4 GHz radio continuum in mjy beam -1 . 


llSauIt et a.1.11 lQQflh . After the initial bandpass calibration us¬ 
ing a WSRT standard calibrator, continuum images are pro¬ 
duced and improved by iteratively applying a continuum 
self-calibration. The resulting calibration tables are then 
used for the line imaging as well. After a second continuum 
subtraction, the resulting Hi image cubes are produced by 
using a robust weighting of 0.4 in order to produce the opti¬ 
mal balance between angular resolution and surface bright¬ 
ness sensitivity. We average across every second channel to 
improve the signal-to-noise sensitivity which results in a ve¬ 
locity resolution of 8.5 kms -1 (after Hanning smoothing). It 
should be noted that these final reduction parameters were 
determined after we had already tested various combinations 
of weighting schemes (including natural and uniform weight¬ 
ings) with different angular and frequency tapering for each 
set of observations. Table[2]lists the RMS levels obtained for 
each galaxy in our sample. 


4 RADIO CONTINUUM PROPERTIES OF 
BLUE EARLY-TYPE GALAXIES 

We observe 1.4 GHz radio continuum emission in the nuclear 
regions of J0900+46, J1117-I-51 and J1237+39. No nuclear 
radio emission was found for J0836+30, but we do observe 
two radio lobes extending 88.4 kpc north-west and 102.5 kpc 
south-east of J0836+30 along the same direction as the ex- 
tragalactic Hi cloud. Figure[5]shows a multicolour composite 
image of the J0836+30 field where the radio continuum and 
the integrated Hi emission is shown in magenta and cyan, 
while the background three-colour image is produced from 
the five SDSS ugriz bands. Further discussion of the pos¬ 


sible mechanism for the displacement of the Hi reservoir in 
J0836+30 can be found in Section 5.1. 

We argue that the two radio lobes flanking both sides 
of J0836+30 are likely to be faded radio lobes belonging to 
J0836+30 and not likely to be due to a chance alignment 
with background sources. Although there is some overlap 
between the South-Eastern radio lobe with a WISE 3.6 mi¬ 
cron source, we do not find an optical/IR counterpart for the 
North-Western radio lobe. In addition, the projected sepa¬ 
ration between both lobes and J0836+30 are similar enough 
for the differences in angular separation (as well as bright¬ 
ness) to be accounted for by projection effects. 

The nuclear radio luminosities range from 2.0 x 10 35 to 
7.4 x 10 35 ergs s _1 (or 1.0 x 10 2 ° to 3.7 x 10 21 W Hz -1 ) and 
are consistent with those found mostl y in other radio studies 
of low-redsh i ft sta r-forming galaxies dMauch fc S adler 2007; 

I Jarvis et al.ll2ofol ). For J0836+30, the luminosities of the 
radio lobes in the north-west and south-east directions are 
9.26x 10 35 and 1.22x 10 37 ergs s _1 , respectively. Table[3]lists 
the measured 1.4 GHz radio continuum properties. FigureQ] 
shows the radio continuum emission as green contours and 
the Hi emission as white contours for the entire pilot sample. 

Assuming that the nuclear radio continuum emission 
is due to star formation, we can estimate the star forma¬ 
tion rate upper limits in our pi lot sample . Using the star 
forma tion rate calibrations from I Belli ([2003 ); Ho pkins et al.l 
d2003h . we find that our galaxies have star formation rates 
that are fewer than 2 solar mass per year. J1117+51 has the 
highest estimated star formation rate with 1.7 solar mass 
per year. 

Archival IRAS 60 /im observations dMoshir fc et al.l 
1 1990h were found for J1117+51 and J1237+39. Using the 
far-infrared to radio correlation (FRC) of low-redshift star- 
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forming galaxies iYun et all (l200lh . we find that the expected 
far-infrared 60 fim luminosity that correlates to our mea¬ 
sured radio luminosity if the radio emission is due solely 
to star formation is log Leo of 9.6 and 9.0 for J1117+51 
and J1237+39, respectively. In comparison to the IRAS log 
Leo values of 9.4 and 9.3, we find that the predicted val¬ 
ues from the FRC are consistent with the idea that the 1.4 
GHz continuum emission originates from star formation for 
both J1117+51 and J1237+39. It should be noted that the 
FRC has a scatter of approximately 0.25 dex and the IRAS 
measurement uncertainties are at the 20% level. 

Although the radio continuum emission observed from 
J1237+39 appears to be extended, it is in fact the result of a 
chance alignment between J1237+39 and a projected back¬ 
ground galaxy, SDSS J123716.92+392921.9, which is 8.5 
times more distant (at redshift z = 0.176) than J1237+39. 
Figure [6] shows that the observed radio emission can be sim¬ 
ply modelled by two point sources. The lack of any regions 
which display a flux excess or a flux deficit in our resid¬ 
ual image (panel c of Figure [6} is consistent with the idea 
that the radio emission from J1237+39 does not have an 
extended diffuse component. 


5 Hi CONTENT OF BLUE EARLY-TYPE 

GALAXIES 

As the four blue ETGs in our sample represent four sepa¬ 
rate stages of evolution (as traced by the UV/optical colours 
and the [OIII]/[H/3] ratios), we observe a correlation with 
their respective Hi gas-to-stellar mass fractions whereby the 
gas fractions decrease as a function of time with increasing 
UV/optical colours and increasing [OIII] / [H/3] ratios. Ta¬ 
ble [3] lists the measured Hi properties. The galaxy at the 
earliest stage of quenching in our sample (i.e. the galaxy with 
the bluest NUV/optical colour and the weakest [OIII]/ [H/3] 
ratio), J1237+39, features a central, undisturbed, rotating 
HI disk. Both spatial and kinematic asymmetries are not 
observed in this galaxy (see Figure]!] and Figure 0. 

In the second bluest galaxy, J1117+51, the Hi gas ap¬ 
pears to be mostly within J1117+51. However the gas mor¬ 
phology is fairly asymmetric and appears offset from the 
optical centre of the galaxy. Assuming the optical velocity 
to be the systemic velocity of each galaxy, we find the Hi gas 
in J1117+51 to be largely blue-shifted by up to 170 kms” 1 
from the galaxy’s systemic velocity (see Figure 0). 

In the redder galaxies and the galaxies with stronger 
[OIII]/[H/3] ratios—indicative of ionisation levels consistent 
with Seyfert activity, the Hi reservoirs of J0836+30 and 
J0900+46 are completely displaced from their respective 
galaxy centres by 14-86 kiloparsecs. In addition, we find 
a significant amount of red-shifted Hi (by « 70 — 80 kms -1 ) 
with respect to the systemic velocity of the individual galaxy 
(see Figure [T}. 

In J0836+30, the entire gas reservoir is offset spatially 
by 1 areminute (projected distance of 30.5 kpc). To approx¬ 
imate the timescale for the displacement of the gas reser¬ 
voirs from J0836+30 and J0900+36, we simply divided the 
projected distances of the gas clouds by the observed Hi 
velocity offsets from the systemic velocities taken from the 
optical observations. For J0836+30, the timescale for the 
gas reservoir to have reached its current projected distance 


is approximately 373 Myr. Similarly, the Hi emission de¬ 
tected near J0900+46 spans between 14 kpc to 86 kpc from 
the optical center of J0900+46 and appears to be redshifted 
by approximately 70 kms^ 1 (Figure [8]o). We estimate the 
displacement timescale to range between 0.2 to 1.2 Gyr. 


5.1 Mechanisms for the removal of Hi 


Support for the hypothesis that the observed Hi clouds in 
the vicinity of J0836+30 and J0900+46 have been removed 
from their parent galaxy come from two main arguments. 
Firstly, the proximity of the gas clouds are very similar to 
the distances between the Milky Way and its high velocity 
clouds (HVCs). Secondly, neither J0836+30 nor J0900+46 
reside in dense environments with close neighbouring galax¬ 
ies which could be responsible for strip ping the gas reser¬ 
voirs from each of these galaxies. In fact. lVerlev et all (|2007l l 
has that found both J0836+30 and J0900+46 to be iso¬ 
lated galaxies using different galaxy isolation metrics. Sev¬ 
eral other previous studies have a lso classified J0836+ 30 to 
be a prototypical isolated galaxy (IVerdes-Monte negro et al.l 
120051 : IStocke et akll2004 lKarachentsevalll973l '). 

I 11 dense galaxy environments, such as clusters or 
compact gal axy groups, interga l actic gas clouds are not 
unco m mon jKilboriiet ah[ ]2000l : lOosterloo fe van Gorkoml 
120051 : iBorthakur et al! 20ldf due to various gravitational 
and hydrodynamic interaction processes suc h as galaxy- 
galaxy inte ractions (Mihos 200-3), harassment (iMoore et al.l 
1 19961. 1998f).ram pre ssure and turbulent viscous stripping 
( Vollmer et al.lf200lf) . 

Ram pressure stripping is unable to shift such large 
amounts of gas using simple Toom re & Toomre pressure 
arguments (e.g. IChung et al] 120071 1. Also, the galaxies in 
this pilot sample do not have nearby galaxy neighbours nor 
do they reside in dense galaxy environment such as clusters 
or compact groups. 

Nearby examples where the majority of the galaxy’s gas 
is displaced via tidal stripping also results in strong mor- 
phological distort io ns (e.g. NGC 4438 in the Virgo Cluster 
llHota et al.ll2007ln . iHota et al] (l2012fl has also identified a 
post-merger system, NGC 3801, where the extraplanar Hi 
resulting from the merger has yet to be encountered by a 
newly-triggered radio jet. As we do not observe any strong 
optical signatures from our galaxies such as stellar tidal tails, 
we can only infer that (l)the tidal features resulting from 
the strong interaction that removed the large gas mass have 
faded; or (2)another physical mechanism is responsible for 
the gas removal. 

In the latter case, we posit that an active central engine 
has the required energy to expel a galaxy’s gas. In the case 
of J0836+30, we observe two radio lobes in the same direc¬ 
tion as the ejected Hi cloud. This is suggestive that at ear¬ 
lier epochs, J0836+30 hosted a radio AGN which may have 
blown out the gas cloud that we observe. This hypothesis is 
consistent with recent studies which have found nuclear ac¬ 
tivity to be more prevalent in star-forming g alaxies relative 
to we akly star-forming or quenched galaxies (iRosario et al] 
I2013H . 

As the displacement timescale for this cloud is esti¬ 
mated to be « 373 Myr, it is conceivable that J0836+30 has 
since evolved to a more Seyfert state (as observed via the 
optical nebular emission). A test of this hypothesis would 
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Q. 


Figure 6. The 1.4 GHz radio continuum emission from J1237+39. Centred on J1237+39, we model the emission as two overlapping 
point sources. The solid blue line in all three panels mark the 3 a flux level of the observed 1.4 GHz emission (see Panel a). Panels b and 
c show the model of two point sources and the residual maps, respectively. 
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Figure 7. Velocity maps and spectra of J1117+51 (panel a) and J1237+39 (panel b). The 3-color velocity maps correspond to the Hi 
emission in the velocity range highlighted in the same color on the corresponding velocity spectrum. It should be noted that the green 
represents the optical velocity which we assume to be the systemic velocity of each of the blue early-types sampled in this paper. 
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Table 3. Radio (1.4 GHz) and Hi properties of blue early-type galaxies. 


Galaxy 

(i) 

Si.4 
(2) 

Lia 

(3) 

Lia 

(4) 

SFR 1.4 

(5) 

■Shi 

( 6 ) 

u Hi 

(7) 

™Hi 

(8) 

M H i 

(9) 

J0836+30 

< 2.5 x 10 -4 

< 3.3 x 10 20 ** 

< 6.6 x 10 34 

< 0.4 

0.078 (0.025) 

7705 

140 

2.0 (0.7) 

J0900+46 

9.7(2.4) x 10 “ 4 

1.5 x 10 21 

3.0 x 10 35 

1.1 

0.410 (0.182) 

8277 

62 

12.3 (5.5) 

J1117+51 

23.4(3.8) x Hr 4 

3.7 x 10 21 

7.4 x 10 35 

1.7 

0.100 ( 0 . 020 ) 

8185 

140 

3.1 (0.6) 

J1237+39 

12 . 0 ( 1 . 6 ) x 10 -4 

1.0 x 10 21 

2.0 x 10 35 

0.7 

2.448 (0.403) 

6096 

245 

40.8 (6.7) 


**Note that the luminosities of the radio lobes in the north-west and south-east directions are 4.63 (±0.44) x 10 21 WHz -1 and 
6.12 (±0.03) x 10 22 WHz -1 , respectively. Col. (1): Galaxy identification. Col. (2): Total 1.4 GHz radio continuum emission in Jansky. 
Col. (3): Radio continuum luminosity (1.4 GHz) in WH z -1 . Col. (4): Radio continuum lu minosity (1.4 GHz) in ergs s —1 . Col. (4): 

Estimated star formation rate using the calibration by I Belli d2003l) ; Hopkins et alj <120031) in Mq year -1 . Col. (5): Integrated Hi 
emission in Jansky. Col. ( 6 ): Hi radial velocity as measured from the mid-point at full-width half-maximum in kms -1 . Col. (7): The 
50% Hi velocity width in kms -1 . Col. ( 8 ): Total Hi mass in units of xlO 8 M 0 . Values in brackets give the estimated errors. 
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Figure 8 . Velocity maps and spectra of J0836+30 (panel a) and J0900+46 (panel b). The 3-color velocity maps correspond to the Hi 
emission in the velocity range highlighted in the same color on the corresponding velocity spectrum. The white solid circle represents 
the position of the blue early-type galaxy relative to the observed Hi emission. It should be noted that the green represents the optical 
velocity which we assume to be the systemic velocity of each of the blue early-types sampled in this paper. 


be the presence of ionized gas adjacent to and between the 
observed Hi cloud and the galaxy. Therefore, in the case of 
J0836+30, we think that the alignment between the radio 
lobes and the extragalactic Hi cloud is more suggestive of 
gas ejection via radio jets than a previous tidal encounter 
whose tidal signatures have since faded in the past 373 Myr. 

Previous observations of AGN-driven gas outflows have 


so far been from studies where the outflows are determined 
from th e kinematic structures of the emitting or absorb¬ 
ing gas llGopal- Krishna fc Irwin 200d: Nesvadba et al.ll2009l : 


iMahonv et al] 2013! ; iMorganti et all l2013j j . However, recent 
observations of 12 CO (1-0) emission from high-redshift radio 
galaxies have found spatial offsets between the molecular 
gas reservoir and their host galaxies on the scales of several 
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tens of kiloparsecs an d aligned in the sam e direction to the 
radio hotspots/lobes jEmonts et al]|2014lh We propose that 
we are witnessing a similar ejection of cold gas in J0836+30. 

In the case of J0900+46, we neither observe strong radio 
continuum emission nor nearby galaxies which could remove 
the gas from this galaxy. 

We posit that the observed extragalactic gas clouds in 
J0836+30 and J0900+46 are outflows and not inflows for 
two reasons. Firstly, if the observed cold clouds were be¬ 
ing accreted, the source of this gas is either a neighbour¬ 
ing galaxy or primodial gas leftover from galaxy formation. 
However, neither galaxies have neighbouring galaxies and 
the accretion of primodial gas at the current epoch is in a 
warm-hot phase rather than t hat of the observed cold Hi. 
Secondly, iBouche et al.l i 2111 2 :■ found that gas outflows are 
usually found at angles nearly orthogonal (> 60 degrees) 
to the plane of the galaxy, whereas inflows have small az¬ 
imuthal angles (< 30 degrees) the galactic plane. Therefore, 
we think that it’s more likely that our observed extragalactic 
Hi clouds originated from our target galaxies. 

The mass of our gas clumps are approximately the to¬ 
tal mass of HVCs around the Milky Way, the LMC and the 
SMC. Assuming that our observed gas clumps are analogous 
to these HVCs, we can expect a similar gas dissipation time 
on the order of several hundred million years if there isn’t 
some sort of support mech anism tha t prevents o r slow down 
against cloud dissipation (jPutman et al.1120121 ). The cloud 
survival time is linked closely to its t otal mass, cloud den - 
sity, relative halo density and velocity dPutman et al.|[2012l) . 
Previous studies have found that while the least bound ma¬ 
terial is likely to expand out into the IGM, bound structures 
are likely to fall back onto the galaxies in less than 1 Gyr 
ijHibbard fc Mihoslll99Rl : iHibbard fc van ( iorkomlll 99(l l. 


5.2 Hi gas fractions 


Using the scaling rela tions fro m the a40 catalog of the AL¬ 
FALFA survey ( Huang et alj 120121 '), star-forming galaxies 
with the stellar masses of J0836+30 and J0900+46 would 
be expected to have average Hi masses of approximately 
1 x 10 10 Mq and 8.7 x 10 9 Mq, yet we observe only 2 x 10 s 
M@ and 1.2 x 10 9 Mq, respectively. This implies that some 
fraction of the gas may have been ionized and heated and 
we are currently only observing the remaining fraction of 
Hi. The external Hi clouds we observe represent the rem¬ 
nants of the original gas reservoirs swept out of the host 
galaxy, and from the present data it is not clear whether the 
missing Hi was also expelled or heated or simply dissipated. 
Regardless, the displaced gas reservoir is required to com¬ 
plete the quenching of these galaxies and prevent further 
star formation. 

Extrapolating from the Kennicutt-Schmidt star forma¬ 
tion law, the correlation between the Hi gas fraction and the 
star formation history traced by the NUV —r col or provides a 
rough indicator for the co ld gas surface density ([Huang et al.1 
120121 ; IZwaan et akll2013li . Relative to gas-rich star-formin g 
galaxies (mostly spirals) (ALFALFA; iHuang et al.1 12012 1. 
low-redshift transition gala xies with ste l lar masses greater 
than 10 10 Mq (GASS; ICatinella et al.l l2012l l and nearby 
early - type galaxies from the Atlas3D survey (Serra et al.l 
120121 : iGappellari et al.1 120131 : lYoung et, al.l 120131 '). we find 
that J1117-51 and ,11237+36 have similar gas fractions and 


NUV —r colors to E+A galaxies and straddle the region be¬ 
tween star-forming galaxies with low gas fractions and that 
of gas-rich early-type galaxies. 

Figure [9] compares the Hl-stellar mass ratio and the 
integrated NUV —r colour of J1117+51 and .11237+39 
(marked as black stars) to those of E+A galaxies (rep¬ 
resented by black crosses) where Hi has been previously 
detected (IChang et al.l l200ll ; iBuvle et al.l 120061 : iHelmboldu 
120071 : IZwaan et al. 20131 '). galaxies from ALFALFA are rep¬ 
resented by blue solid points and the solid line marks the av¬ 
erage colours and average gas fractions found from GASS. 
Near-ultraviolet (NUV) measurements of our sample were 
obtained from the Galex satellite telescope via the NASA 
Extragalactic Database^. 

We also compare our pilot sample of blue ETGs to re¬ 
cent Hi observations o f nearby ETGs from the Atlas3D sur¬ 
vey (iSerra et al.ll2012i ). The Atlas3D ETGs with regularly- 
rotating undisturbed Hi morphologies are represented by 
open circles, while, Atlas3D galaxies with very disturbed 
and unsettled Hi morphologies and kinematics are repre¬ 
sented by solid diamonds in Figure [9] It should be noted 
that all the Atlas3D ETGs with disturbed Hi morphologies 
reside in group or Virgo Cluster environments. Also, for a 
given NUV —r colour, these Atlas3D galaxies have higher 
gas fractions than that of the E+A or blue ETG samples. 
This suggests that either (1) the Atlas3D ETGs where Hi is 
observed are at earlier stages of evolution than the sample 
of E+A galaxies or blue ETGs; or that (2) the merger sce¬ 
nario is as likely to result in an increase in gas fraction as it 
is to a reduction. Upper limits are shown for J0836+30 and 
J0900+46 due to the lack of Hi within the host galaxies. 

We find the Hi gas fraction and NUV —r colour of 
J1237+39 (the galaxy at the earliest stage of quenching) to 
be very comparable to those of star-forming galaxies from 
both the ALFALFA and GASS surveys. On the other hand, 
the gas fractions for the three more evolved blue ETGs 
in our sample are significantly lower than the average gas 
fractions expected from star-forming galaxies, transitioning 
galaxies or even gas-rich early-type galaxies with similar 
NUV —r colours. Similar to our results, the gas fractions of 
post-starburst E+A galaxies are lower than the average gas 
fractions for any given NUV —r colours. The combination of 
the Hi mapping and the radio continuum observations of 
J0836+30 and J0900+46 suggest that the star formation in 
these two galaxies wi ll be truncated so o n. Th is result is con¬ 
sistent with those of ISchawinski et al.l (l2009l l who were un¬ 
able to detect any significant molecular gas reservoirs within 
Green Valley galaxies with Seyfert ionisation properties. 

As previously seen in Figure [3] these pilot Hi observa¬ 
tions show that the main mechanism for the fast quenching 
of star formation in blue early-types is due to the physical 
displacement of the main gas reservoir from which stars are 
formed. Hence, it is likely that the depressed gas fractions 
from our sample (see Figure [9]) are due to a fast quenching 
process (relative to other evolving galaxies) which removes 
a significant fraction, if not the entire gas reservoir, from 
which stars are formed. 

Should the gas fall back, it may restart minor star 
formation which may be visible as the frosting « 1% 


2 http://ned.ipac.caltech.edu 
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Figure 9. Hi-to-stellar mass ratio as a function of galaxy color. 
The solid line shows the average gas fracti ons found from GAS S 
survey of massive transition-type galaxies (|Catinella et al.ll2012h . 


mass fraction of yo ung stars ofte n observed in quenched 
early-type g al axies (| Yi et al.l 1200, r l ISchawinski et all l200fil : 
iKavirai et al.l 120071 ). While it may cause frosting, this re¬ 
turning gas however will not return the host galaxy to self- 
regulated star-formation on the main sequence. An alter¬ 
native fate for the ejected gas is that it persists at large 
distances from the quenched galaxy, as is observed in the 
Milky Way (in the form of the starless Magellanic stream; 
iFor et all 2014h : as well as in many quenched early-type 
galaxies fe.g. fSerra et ah ISmT 


5.3 Argument against coincident dark galaxies 

Dark galaxies can be typically defined as the extreme end 
of low surface brightness galaxies where few stars are found, 
consisting mainly of gas. With respect to our observed extra- 
planar gas clumps, it is very unlikely that these gas clumps 
are neighbouring dark galaxies. Current re-ionization mod¬ 
els of the Universe predict the latter as they find that the 
gas from 95 per cent of the low-mass systems (M v iriai ^ 10 s 
Mg or Voire 20 km s -1 ) appears to have been photoevap- 
orate d during the epoch of re-ionization (ISusa fe Umemural 
|2 oo4i . 

Apart from a few reported cases of extragalactic 
gas clouds in group/clu ster environments dMinchin et al.l 
120051 : IDavies et all 12006m with no known optical counter¬ 
parts, c urrent and previous Hi all-sky su rveys such as AL¬ 
FALFA (lllavne s et al . 20nh and HIPASS (IWong et al.ll2009l: 
iDovle et al.l 120051 1 have not found any isolated extragalac¬ 
tic HI clouds devoid of stars. This is consistent with simple 
gas equilibrium models ( iTavlor fe Websteill2005l j which con¬ 
cluded that in the absence of an internal radiation field, dark 
galaxies/gas clouds with masses greater than 10 9 Mg will 
become Toomre unstable against star formation and start 
forming stars. Galaxy evolution simulations demonstrated 


that the majority of gas clouds identified with no known op¬ 
tical counter part can be easily reproduced by galaxy-g alaxy 
interactions llBekki et alJl2005l : iDuc &; Bournaudll2008l l. 


6 SUMMARY 


We have performed deep imaging of the Hi content and 1.4 
GHz radio continuum emission of four blue early-type galax¬ 
ies that are at four different stages of star formation trunca¬ 
tion using the WSRT. The Hi morphologies, kinematics and 
relative Hl-gas fractions are excellent probes and measures 
of the quenching evolutionary stages for each of our sample 
galaxy. A summary of our results are as follows: 


(i) We observe nuclear 1.4 GHz radio continuum emission 
that are consistent with emission from star formation from 
three galaxies at the three earliest stages of quenching evo¬ 
lution, namely, J1237+39, J1117+51 and J0900+46 (in the 
order of earliest to later stages of evolution). 

(ii) The galaxy at the earliest stage of quenching, 
J1237+39, also has the bluest NUV —r colour, weakest 
[OIII/H/3] ratio, the highest Hi gas fraction and a symmetric 
rotating Hi disk. 

(iii) At more advanced stages of quenching evolution, we 
observe increasingly asymmetric and increased spatial off¬ 
sets between the Hi gas and the stellar component of the 
galaxy. Non-rotating Hi gas kinematics are also observed. 
In the case of the two galaxies at the most advanced stages 
of evolution where Seyfert ionisation signatures have been 
observed (J0900+46 and ,10836+30), the Hi gas reservoirs 
have been entirely expelled by approximately 14-86 kpc 
from their respective host galaxies. Due to the lack of neigh¬ 
bouring galaxies, it is difficult to attribute the stripped gas 
to tidal or ram pressure interactions. 

(iv) In the galaxy at the most advanced quenching stage, 
J0836+30, the expelled gas reservoir is observed to be in 
alignment between the host galaxy and two radio lobes— 
suggesting the gas reservoir may have been swept out of 
the galaxy by powerful outflows from the central AGN in a 
previous active phase. This scenario is consistent with recent 
observations of AG N-driven gas outflows both locally and 
at higher redshifts ([Sturm et al.l 20 111 : Mahonv et al.ll2013l : 
iMorganti et al.l 12014 Emonts et all 20141') . 

(v) We conclude that the rapid quenching of star forma¬ 
tion in low-redshift early-type galaxies is due to the physical 
expulsion of the entire gas reservoir, rather than the exhaus¬ 
tion of gas via star formation. An AGN outflow has the nec¬ 
essary energy to expel the reservoir of cold gas. In addition, 
we do not think that AGN heating of the gas is a dominant 
mechanism because the Hi gas would not remain visible as 
a coherent structure if this was the case. 
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